Cosmic Star Formation History to z=l from a Narrow Emission 
Line Selected Tunable Filter Survey 
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ABSTRACT 

We report the results of a deep 3D imaging survey of the Hubble Deep Field North using the 
Taurus Tunable Filter and the William Herschel Telescope. This survey was designed to search 
for new line emitting populations of objects missed by other techniques and to measure the cosmic 
star-formation rate density from a line-selected survey. We observed in three contiguous sequences 
of narrow band slices in the 7100, 8100 and 9100A regions of the spectrum, corresponding to a 
cosmological volume of up to 1000 Mpc 3 at z = 1, down to a flux limit of ~ 2 x 10 -17 ergs 
cm~ 2 s^ 1 . The survey is deep enough to be highly complete for low line luminosity galaxies. 
Cross-matching with existing spectroscopy in the field results in a small line-luminosity limited 
sample, with very highly redshift identification completeness containing seven [Oil] , H/3 and Ha 
emitters over the redshift range 0.3 - 0.9. Treating this as a direct star-formation rate selected 
sample we estimate the star- formation history of the Universe to z = 1. We find no evidence for 
any new population of line emitting objects contributing significantly to the cosmological star- 
formation rate density. Rather from our complete narrow-band sample we find the star-formation 
history is consistent with earlier estimates from broad-band imaging surveys and other less deep 
line-selected surveys. 

Subject headings: Galaxies: star formation, Galaxies: evolution 



1. Introduction 

Studies of cosmic star formation history over 
the past decade have been an integral part of the 
development of our knowledge of galaxy forma- 
tion and evolution. These studies have found that 
star formation rates rise significantly from z = 
to z = 1, peak between 1 < z < 2, and flat- 
ten or decline (depending on the role of dust ob- 
scuration on ultraviolet measures) beyond z = 2 
(Madau et al. 1998; Connolly et al. 1997; Stei- 
del et al. 1999; Glazebrook et al. 1999; Flores et 
al. 1999; Giavalisco et al. 2004). Star formation 
rates (SFR) are usually determined by measuring 
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the luminosity densities of galaxies, either from 
emission lines or from the continuum spectrum of 
a wavelength range. Highly luminous stars live 
very short lifetimes, and can be used as a proxy 
for instantaneous star formation rates in galax- 
ies. These stars easily ionize the gases surrounding 
them, and the strength of the recombination lines 
of various elements gives the number of ionizing 
photons, which can be converted into star forma- 
tion rates (Kennicutt 1992). They also give rise to 
broad-band ultraviolet emission which can also be 
used as a measure of star-formation rate accessi- 
ble at high-redshift (Madau et al. 1998), but this 
should be used with caution as the ultraviolet is 
much more subject to dust extinction than the op- 
tical re-combination lines and high-redshift SFRs 
can be seriously underestimated (Glazebrook et al. 
1999; Steidel et al. 1999). 

In the intermediate redshift < z < 1 regime 
the most commonly studied emission lines for mea- 
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suring the cosmological SFR are Ha (Gallego et al. 
1995; Tresse and Maddox 1998; Glazebrook et al. 
1999; Tresse et al. 2002; Kewley et al. 2002; Fujita 
et al. 2003), and the 3727A Oil line (Hogg et al. 
1998; Kewley, Geller, & Janscn 2004). Ha, when 
corrected for extinction is directly proportional to 
the number of ionizing photons from massive stars 
(Osterbrock 1989). Since such stars are short-lived 
(~ 20 Myr) their abundance traces the instan- 
taneous SFR. H/3 can also be used this way but 
the extinction and stellar absorption corrections 
are larger. The [Oil] emission line (3727A) has a 
more complex dependence on metallicity (Kewley, 
Geller, & Jansen 2004) and is less straight-forward 
to use as a SFR tracer. 

Most studies use broad band filters to select 
a magnitude and/or color- limited samples which 
gives rise to selection effect issues when computing 
the total SFR. For example populations of strong- 
lined high equivalent width objects may contribute 
significantly to the SFR while having weak con- 
tinuum emission and hence lying below typical 
broad-band magnitude limits. An alternative ap- 
proach is to use a narrow-band filter to directly 
select emission line objects down to an emission 
line flux limit. So far as the line luminosity is in 
proportion to the SFR then this is equivalent to 
a SFR-selected sample and if the line's redshifted 
wavelength matches that of the narrow filter's one 
obtains a higher signaknoise on lower luminosity 
objects in a given integration time. Also, this tech- 
nique produces a volume-limited sample, since the 
narrow observed bands correspond to small win- 
dows in rcdshift space with approximately con- 
stant luminosity limits. 

Previous narrow-band filter surveys have al- 
most universally relied on using special interfer- 
ence filters. Some of the deepest, most recent work 
includes: Pascual et al. (2001) who surveyed 684 
arcmin 2 to a flux limit of about 5 x 10~ 16 erg s _1 
cm~ 2 , finding 16 objects at z = 0.24 and Moor- 
wood et al. (2000) surveyed 100 arcmin 2 in the 
near-infrared to a flux limit of 5-12 x 10 -17 erg s _1 
cm~ 2 , finding 6 objects at z — 2.2. A recent target 
of these searches has been high-z Lyman-a emis- 
sion: Steidel et al. (2000) surveyed 94 arcmin 2 in 
a deep Keck exposure reaching fluxes of 4xl0~ 17 
erg s _1 cm~ 2 ; Hu et al. (2004) extending earlier 
work (Hu, Cowie, & McMahon 1998; Hu, McMa- 
hon, & Cowie 1999) surveyed 0.25 deg 2 to a depth 



of 2x10" 17 erg s" 1 cm" 2 and Rhoads et al. (2000) 
reached the same flux limit over 0.72 deg 2 ; how- 
ever spectroscopic follow-up of these surveys has 
been limited to z > 3 objects (winnowed by cither 
broad-band color and/or high-equivalent width se- 
lection). 

A novel approach to selecting line emitters at 
high-redshift is to use a tunable filter, of which 
the most practical technology to realise this is 
based on the Fabry- Perot etalon (Bland-Hawthorn 
2000). Tunable filters have unique advantages 
starting with the fact that a contiguous scan can 
be made in wavelength/redshift space allowing a 
line to be unambiguously identified via a peak in 
flux without reference to any broad-band image. 
They also offer much higher resolution (AA ~ 
10 A) than interference filters (AA ~ 100 A) which 
better match the typical line widths of galax- 
ies with consequent signal: noise gains. The scan 
through the line also permits any continuum to be 
precisely removed in the analysis allowing accurate 
integration through the line profile to measure the 
flux. This flux can also be integrated over the en- 
tire galaxy (essentially 3D aperture photometry) 
meaning total line fluxes can easily be measured 
without the dubious aperture corrections required 
for narrow slit spectra. 

There have been few surveys using this tech- 
nique. Jones and Bland-Hawthorn (2001) using 
the Taurus Tunable Filter on the Anglo- Australian 
Telescope surveyed 972 arcmin 2 to a flux limit of 
5-10 x 10~ 17 erg s _1 cm~ 2 , finding 696 objects 
in six bands from z — 0.08 to z — 0.4. Hippelein 
et al. (2003) surveyed 400 arcmin 2 to a flux limit 
of 3 x 10~ 17 erg s _1 cm~ 2 , finding 438 objects 
at 0.25 < z < 1.4. Neither survey has extensive 
spectroscopic follow-up to give unambiguous iden- 
tification. Jones and Bland-Hawthorn (2001) ar- 
gued based on model luminosity functions that at 
their bright flux limits Ha emitters would domi- 
nate the counts; there was a small amount of over- 
lap (18 of the brighter objects) with the Autofib 
redshift survey (Ellis et al. 1996) which backed up 
their conclusion. Hippelein et al. (2003) had ad- 
ditional observations in a set of 15 medium and 
broad band filters which allowed the redshifts to 
be constrained using fits to the spectral energy dis- 
tributions. Spectroscopic followup of ~ 50 objects 
with R < 25 showed this was reliable for ~ 80% of 
galaxies; other emission line objects without con- 
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tinuum detections amount to < 7% of this sample 
(Hippelein 2004, private communication). 

In this work we report on a small area (20 
arcmin 2 ) but extremely deep (1.7 - 2.4 x 10~ 17 
erg s -1 cm~ 2 ) survey carried out with the Tau- 
rus Tunable Filter on the William Hcrschcl Tele- 
scope in the Hubble Deep Field North (HDF-N) 
(Williams et al. 1996). This field has the advan- 
tage of extensive spectroscopy (Cohen et al. 1996, 
2000); so while the resulting sample is small it is 
unique in being almost completely spectroscopi- 
cally identified. The goal was to explore the pos- 
sibility that additional star formation, which could 
contribute strongly to the cosmological SFR den- 
sity, could be occurring in galaxies with strong 
emission lines but little or no continuum emission. 

The plan of this paper is as follows: in Sec- 
tion 2 we review the details of the observation, 
data reduction, object identification and catalog 
making. In Section 3 we discuss our process for 
calculating galaxies' SFRs and correcting for ex- 
tinction. In Section 4 we describe the method- 
ology we followed for determining the luminosity 
density based on the narrow band line fluxes we 
observed, and calculate the corresponding star for- 
mation rates. In Section 5 we present a cosmic 
star formation history, calculated from our sample 
to z=l, and discuss the implications. Through- 
out this paper we use a cosmology: Ho = 70 km 
sec -1 Mpc" 1 , fl m =0.3, and Q\ —0.7 motivated by 
Spergel et al. (2003). 

2. Observations and Data Reduction 

The observations were carried out on 1997 
March 11-14 at the 4.2m William Herschel Tele- 
scope (WHT) in the Canary Islands using the 
Taurus Tunable Filter (Bland-Hawthorn & Jones 
1998) on loan from the Anglo- Australian Obser- 
vatory. The TAURUS2 instrument (Unger et al. 
1990) on the WHT was identical to the one on 
the Anglo- Australian Telescope. 1 Operating at 
//4 the instrument delivered a pixel scale of 0.29 
arcsecs with the 1024x1024 Tektronix CCD. Con- 
ditions were photometric for most of the run with 
seeing of ~ 0.7arcsec. 

The field centre was RA 12 h 36 m 46.8 s DEC 



Both TAURUS2 instruments have now been decommis- 
sioned. 



+62° 14' 46" (J2000) and the WHT images in- 
cluding the whole of the Hubble Deep Field North 
together with its surrounding environs. By choos- 
ing a TTF etalon gap and broad-band blocking 
filter one generates a slice in wavelength (which 
we call a 'channel') to observe. We scanned the 
etalon in 3 contiguous wavelength regions, each 
a subset of the wavelength range defined by 3 
blocking filters: II (7070 A center/260 A wide), 15 
(8140/330A) and 18 (9090/400A). These are cho- 
sen to lie in the darkest regions of the night sky 
spectrum between 7000A and lOOOOA. The wave- 
length step for the scan was chosen to match, 
approximately, the Full- Width Half-Max of the 
etalon resolution. For etalon Lorentzian profiles 
this provides adequate sampling. We scanned (i.e. 
changing the etalon gap to produce different filter 
wavelengths) through 10 channels in the II and 18 
filters and through 5 channels in the 15 filter. The 
full wavelength scanning parameters together with 
the exposure times per position are given in Ta- 
ble 1. Longer exposure times were used in the red- 
der filters, which correspond to higher redshifts for 
a given line, to partially counteract cosmological 
dimming, The deepest scans were in the 18 filter 
where we exposed for 1.5 hours in each channel. 

2.1. Data Reduction 

The reduction and analysis of tunable filter im- 
ages is discussed in detail in Jones, Shopbell, & 
Bland-Hawthorn (2002). The goal in our case was 
to produce flux calibrated, background-subtracted 
data cubes (i.e. RA, DEC, AA). The raw CCD 
frames we de-biased and flat-fielded using dome 
flats. Frames taken at the same wavelength were 
stacked with a cosmic-ray rejection filter. The 
overall wavelength calibration of the etalon gap 
was determined by taking a fast scan of an Ar- 
gon arc lamp using just the central part of the 
CCD windowed. In the full frame data images the 
Jacquinot spot is approximately centered on the 
CCD; the etalon phase shifts as one goes off-axis 
causing the tuned wavelength for a given chan- 
nel to vary bluewards up to an additional 10 A. 
This effect also gives rise to a broad ring-shaped 
structures in the CCD background; the rings cor- 
respond to the position of night-sky lines in the 
blocking filter band-pass. This was removed by 
fitting for the centre of the ring pattern and then 
subtracting the median count in annular bins. 
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The remaining background after this step still has 
structure clue to CCD fringing; this is more ran- 
dom but had structure on ~ 20-30 pixel scales. 
We attempted to remove this by subtracting a lo- 
cal median, in a matching box, from each pixel. 
This procedure worked, essentially perfectly, in 
the II and 15 filters but was imperfect in the 18 
filter due to the much greater fringing at 9100A. 
This gives rise extra non-Poisson noise which must 
be considered carefully when doing object detec- 
tion (see Section 2.2). 

The overall flux calibration was established by 
taking a fast windowed scan of a spectroscopic 
standard star and a bright reference star in the 
HDF-N field during photometric conditions. This 
was then referenced to the counts of the same ref- 
erence star in the final data cube (which included 
some slightly non-photometric data) in order to 
establish the final calibration of counts into flux 
versus wavelength. 

2.2. Object detection and completeness 
tests 

Our goal was to construct a true three-dimensional 
catalog so we did not wish to require that an object 
have detectable flux in every wavelength channel. 
Rather if an object was detected in one chan- 
nel we wanted to have reliable flux limits in all 
other channels based on the aperture defined by 
the detection channel. To achieve this we ran the 
SEXTRACTOR object detection software (Bcrtin 
& Arnouts 1996) on each 2D channel separately 
using a threshold of at least 5 pixels > 1.2 a above 
the background using pre-smoothing with a Gaus- 
sian kernel matched to the seeing. The noise a 
was set by the CCD readnoise and gain param- 
eters. Certain parts of the 2D images (~ 17%) 
had residual contamination features due to opti- 
cal ghosts and bad columns, these were masked 
off by a binary SEXTRACTOR mask to avoid 
object detection in these areas. The effective area 
searched for objects was then 20.35 arcmin 2 . 

From each individual channel catalog objects 
were photometered in all other channels using the 
same SEXTRACTOR aperture parameters (i.e. 
'two-image' mode for all other channels as the sec- 
ond image). This process was then repeated for 
each channel. What we end up with is a large set 
of catalogs, one for each channel, containing the 
spectra (i.e. aperture flux versus channel) of each 



object based on it's detection in each channel. The 
next step was to merge these catalogs to remove 
duplicates where bright objects were detected in 
many channels. To detect duplicates a centroid 
proximity test of < 2 pixels was used. To remove 
duplicates we kept only the spectrum based on the 
channel with the strongest detection. Some addi- 
tional cosmic-ray filtering was performed at this 
step; essentially removing any source too sharp to 
be a real object. The result of the merging proce- 
dure was a master catalog of all possible objects 
(351 total) in all possible channels with associated 
spectra and noise vectors. This 3D approach is 
very similar to that of Jones and Bland-Hawthorn 
(2001) but differs from Hippelein et al. (2003) who 
used the sum of the channel images for object de- 
tection. For objects which only appear as a line in 
one channel the latter approach is less sensitive. 

The final step is to search for emission lines in 
the final spectra. Choosing the right signaknoise 
cut is critical to avoid contamination by spurious 
objects. One test we performed was to compare 
the actual variance in the spectra with that calcu- 
lated based on CCD gain and readnoise. (The ma- 
jority of objects will not have a spectral feature in 
the observed wavelength range as it is so short and 
so the variance is dominated by noise). We found 
these agreed well in the II and 15 channels but fell 
short by a factor of 2 in the 18 channel. This was 
because 18 had additional channel-to-channel fluc- 
tuations due to significant residual fringing. To 
compensate for this we forced a match by dou- 
bling the calculated noise in the 18 channel. We 
estimate the final 1-er noise in each channel to be 
(7.9,6.0,5.5) x 10~ 18 ergs cm" 2 s" 1 for II, 15 and 
18 respectively. 

For emission line detection we smoothed the 
spectra with a 3 pixel kernel (values [0.33, 1, 0.33]) 
in order to approximately match the Lorcntz 
profile of unresolved lines). The continuum, if 
present, was defined by the median flux value in 
each filter. A line was 'detected' if any pixel had 
flux > 4.5cr above the continuum, this threshold 
being equivalent to 3cr in the unsmoothed data. 
This gave a '3<t' emission line catalog of 41 'ob- 
jects'. Inspection of the images showed that many 
of these objects looked more like artifacts than 
real objects; the effect of the fringing being to im- 
part a non-Gaussian characteristic to the noise.We 
decided to calibrate the 'believability threshold' 
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Fig. 1 . — Narrow band images and resulting spectra of our sample of 8 emission line galaxies (top - bottom, 
order matching Table 2). The left panel shows the individual channel images (labeled with the channel 
number) in the vicinity of the detected line, with the center column corresponding to the peak flux where 
the line is detected. The peaking of the flux is obvious, even in the 2D images. The right panel shows the 
full channel spectra, i.e. flux vs channel number (0-24). The arrow denotes the emission line location in 
each spectrum. The thick solid vertical lines represents the boundary between different filters where the 
wavelength scale is discontiguous. The thick solid horizontal lines denote the measured continuum level in 
each filter. Filter labels are given. 



empirically by searching for emission line in 1000 
randomized sky positions (constrained to be dis- 
tant from all known objects in the field). This let 
us estimate how often by pure chance, given the 
skew noise distribution in our images, we would 
incorrectly detect an emission line when photome- 
tering our data. Spectra were extracted in 3 arcsec 
diameter apertures and the same detection algo- 
rithm applied to measure the false-positive rate. 
At 3<T the contamination due to non-Gaussian 
artifacts is severe: in our master catalog we cal- 
culated a rate of 40 false positives for every 351 
random positions in our synthetic catalog. Since 
the actual catalog of 351 objects had 41 emission 



line detections this led us to conclude that at this 
S/N threshold we are overwhelmingly dominated 
by false objects. Using these simulations as a 
guide we raised the threshold to make a '5er' cat- 
alog of 8 objects — the random catalog showed in 
this we expect on average only 0.7 objects to be 
false. This rapid drop-off in the number of 'ob- 
jects' in the random catalog when moving from 3<r 
to 5cr (only a factor of 1.66 in flux) is consistent 
with the objects in the former catalog being noise 
artefacts and not serendipitous galaxies. Visual 
inspection confirmed that these 8 objects were 
very likely genuine galaxy images. We show im- 
ages and spectra of all 8 objects in Figure 1 and 
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tabulate them in Table 2. 

2.3. Redshift identification 

Rcdshifts were confirmed by cross-matching the 
sky-coordinates with published HDF spectroscopy 
from Cohen et al. (1996, 2000). In 7 out of the 8 
cases there was a clear, positive, unique identifi- 
cation where a known emission line galaxy at the 
correct RA,DEC had the correct redshift to put a 
strong line in the right filter. (The probability of 
the last happening by chance is of order 1%). Out 
of these 7 objects one is in the HDF itself (#8) 
and the other 6 arc in the flanking fields. The 8th 
object is also in the flanking fields and is blended 
with a neighboring object in ground-based images 
of Capak et al. (2004). Its identification remains 
a mystery but we note the object is detected in 
the B-band so it is unlikely to be z — 5.5 Lya. 
We further note the peak flux is at the red end of 
the 18 scan so it could correspond to a continuum 
break instead of an emission line. 

The final rcdshifts and line identifications are 
given in Table 2. Figure 1 shows both the data 
cube images of each object in the vicinity of the 
line and resulting channel spectra. The emission 
lines are clearly evident in both. Continuum was 
also significantly detected in these objects and the 
flux levels agree with the published /-band HDF-N 
photometry of the same objects. The magnitude 
range of the objects is 18.0 < I < 20.2. Since 
the number of objects was small and the magni- 
tudes were relatively bright our main result — no 
surprises — is presaged here. However the line 
flux limits reached are quite faint (see above and 
Table 3): in these redshift windows probed this 
sample represents all the objects down to a small 
fraction of the Schechtcr luminosity L*. 

3. Line luminosity determination and Ex- 
tinction Corrections 

Our method to calculate SFRs is to calculate 
from each observed line luminosity the equivalent, 
extinction-corrected, Ha luminosity and then con- 
vert Ha to SFR using a constant converstion ratio 
which we assume is independent of redshift (this 
is very similar to assuming a constant IMF). 

To standardize our extinction correction across 
the three different observed lines, we converted the 
H/3 and [Oil] fluxes to equivalent Ha fluxes, then 



corrected the Ha fluxes for extinction. We have no 
alternative but to assume a mean extinction in or- 
der to make progress. To dust-correct the H/3 flux, 
we use the relation derived by Tresse and Maddox 
(1998) which covers a similar redshift range to our 
H/3 and Ha galaxies: 



~* Hf3,raw 



2.86 x 10" 



-C(-0.323) 



(1) 



where we use their mean C = 0.46, which corre- 
spondes to A v = CR/1A7 = 1.0 mag assuming 
R=3.2 (Scaton 1979). 

For [Oil], we use the observed ratio: 



Hcx.raw 



FoiI,r 



= 2.22 



(2) 



from Kennicutt (1992). There may be more com- 
plex dependencies of this ratio on, for example, 
luminosity or mctallicity (e.g. Jansen, Franx, 
& Fabricant (2001)). In particular more recent 
analyses of much larger local samples (Jansen, 
Franx, & Fabricant 2001; Hopkins et al. 2003) have 
found F Ha ^ raw / F n. raw = 4.34). Since we have 
one [Oil] object in our sample it is only neces- 
sary to have one conversion factor, we choose to 
use the Kennicutt (1992) value. This also agrees 
with the z = 1 value found by Glazcbrook et 
al. (1999); however that work concerned luminous 
L* galaxies whereas our object is somewhat sub- 
L*. If a higher ratio was adopted this would in- 
crease our highest redshift measurement of the 
star- format ion rate density (below) by a factor of 
two. Our values are currently at the low-end of ob- 
served values and so this would bring the results 
more in to agreement. However given the large 
error bar from one object this point is somewhat 
moot. 

To correct the equivalent Ha fluxes for extinc- 
tion, we adopt the Ay = 1.0 mag value giving 



Fhci, corrected — 2.12 X Fu a . 



(3) 



The raw and corrected equivalent Ha luminosi- 
ties are shown in Table 2. We then use the Kenni- 
cutt (1998) conversion SFR = L(Ha) x 7.9 x 10~ 42 
Mq yr -1 / ergs cm 2 s _1 (which assumes a Salpeter 
(1955) IMF) to convert to a star formation rate. 

Figure 2 shows how the resulting line lumi- 
nosities we measure for individual galaxies com- 
pare with others in the literature ([Oil] and Ha 
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Table 1: Etalon scan parameters in the three different filters. Each wavelength range and emission line 
identification defines a redshift range and corresponding volume, shown below. Volumes are in units of 
Mpc 3 . 



Parameters of 
Etalon scan 


Exp/channel 
(Sees) 


Ha 

z 


Volume 


H/3 

z 


Volume 


[OH] 

z 


Volume 


7058-7112A, 10 steps of 
AA = 6.0A, FWHM = 6.3A 


800 


0.075-0.084 


7 


0.452-0.463 


197 


0.894-0.908 


599 


8100-8133A, 5 steps 
of AA = 8.3A, FWHM = 8.9A 


1800 


0.234-0.239 


31 


0.666-0.673 


203 


1.173-1.182 


462 


9041-9125A, 10 steps of 
AA = 9.3A, FWHM = 13.5A 


5400 


0.378-0.390 


176 


0.860-0.877 


696 


1.426-1.449 


1340 



luminosities from Tresse et al. (2002), Hippelein 
et al. (2003), which we also standardize to dust- 
corrected Ha luminosities in the same way as for 
our data). Our objects are all only moderately 
star-forming, typically only a few M Q yr _1 and 
our detection limits correspond to quite low lumi- 
nosity values. 

4. Calculation of luminosity and SFR den- 
sities 

To calculate SFR densities (i.e. SFR per unit 
comoving cosmological volume) we use the same 
SFR/luminosity conversion procedures but this 
time applied to the dust-corrected Ha luminosity 
density. Because our flux limits are well below L* 
at all redshifts the very simplest estimator for the 
luminosity density pv- 

M*fc)=£^r (4) 

i 

where 14 is the volume in slice k given in Table 1 
comes very close to estimating the total luminosity 
density irrespective of the shape of the luminosity 
function at the faint end. (Note we are using dust- 
corrected equivalent Ha luminosities). Any new 
population of faint line emitting galaxies would 
show up as an excess luminosity density compared 
to other surveys. 

In order to estimate errors and the amount of 
missing light below our flux limits we next assume 
the Ha luminosity function takes the form of a 
Schechter (1976) function 

<j>{L)dL = <j>*(L/L*) a exp(-L/L*)d(-L/L*) (5) 



and estimate (f>*,L* as follows. Writing x = 
Lum(zk) I L* where Lu m (zk) is the luminosity limit 
corresponding to the flux limit in a given slice, the 
number of galaxies per unit volume n we expect 
to see is: 

/•OO 

n = <j>* x a e- x dx = <t>*T{a + l,x) (6) 

J X 

where T is the incomplete Gamma function. Sim- 
ilarly the integrated luminosity density pl above 
the flux limit is 



p L = (t>*L* / x a+1 e- x dx = (t>*L*T(a + 2, x) 

J X 

(7) 

Then we can estimate L* using: 

T(q + l,x) = n_ 
L*T(a + 2,x) p L [ ' 

which can be solved iteratively for L* given the 
observed values of n and in a given redshift 
slice (in effect wc constrain <j>* and L* to give the 
observed luminosity density > our flux limit, but 
they do play a role in the error calculation below). 
Then <j)* can be derived using equation 6 and the 
total luminosity density calculated by putting x = 
in equation 7. 

The amount of light missed below our flux lim- 
its depends mainly on a. For our final, 'total', val- 
ues we adopt a of —1.3 which is the mean of local 
optical-IR surveys (Gallego et al. 1995; Blanton 
et al. 2001; Huang, Glazebrook, Cowie, & Tin- 
ney 2003). With this value we find that we are 
seeing between 89% of the light in the z = 0.38 
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Table 2: Raw and extinction-corrected Ha equivalent luminosities of each galaxy. Luminosities are shown as 
log L in units of erg s _1 . 



Obj# 


z 


Line 


Raw Line 
Luminosity 


Corrected 
Ha Luminosity 


RA 


Dec 


45 


0.378 


Ha 


40.94 


41.27 


12 37 04.65 


62 16 52.1 


190 


0.378 


Ha 


40.53 


40.85 


12 37 04.04 


62 15 23.2 


60 


0.378 


Ha 


40.83 


41.16 


12 36 39.72 


62 15 26.2 


8 


0.454 


H/3 


40.65 


41.58 


12 36 42.93 


62 12 16.3 


33 


0.457 


H/3 


41.02 


41.95 


12 36 58.39 


62 15 48.7 


47 


0.455 


H/3 


40.71 


41.64 


12 36 31.16 


62 12 36.2 


88 


0.904 


[Oil] 


41.41 


42.08 


12 36 38.80 


62 15 47.2 


330 


? 


? 


? 


? 


12 36 43.18 


62 16 25.1 



bin to 61% of the light in the highest z = 0.90 
bin. If we a were higher we would of course miss 
more. For example if we used a = —1.6, which 
is seen in the ultraviolet luminosity function of 
z > 3 star-forming Lyman break galaxies (Stei- 
del et al. 1999), then our corresponding luminosity 
completeness is only 76-43%. 

Errors are tricky to calculate in samples with 
small numbers of objects (Gehrels 1986), we adopt 
a Monte-Carlo approach where we populate a 
model universe according to the Schechter func- 
tion (equation 5). This lets us account for scat- 
ter both in the luminosities and the numbers of 
objects in a finite survey. A grid of (<f>*,L*) val- 
ues are explored bracketing our measured values, 
for each grid point we make 1000 realizations of 
our observed volume. For each realization we re- 
peat our estimation of the luminosity and number 
densities using Equations 6-8. We can then com- 
pute a mean (p, n) and variance (Cp,cr^) for these 
quantities and form a x 2 value to characterize the 
goodness of fit of the grid point to the observed 
data (p, n): 



x 2 = 



P 



-\ 2 
P 



+ 



(9) 



Then the range of luminosity densities corre- 
sponding to the Ax 2 = 1 contour gives us our 
error bars. Note these are not symmetric and the 
direction of the asymmetry changes with redshift. 
This is because our x 2 statistic has both luminos- 
ity and number components whose relative impor- 
tance changes; this is compounded with the high 



asymmetry of the Schechter function. We note 
that for the z = 0.90 bin where there is only one 
object the uncertainty is a factor of two which is 
in accord with intuitive expectations. The final 
luminosity densities along with our calculated </>* 
and L* are listed in Table 3. We also tabulate the 
results of the simple direct estimate (equation 4) 
as well as our more sophisticated method. 

The same Kennicutt (1998) conversion factor 
from in Section 3 converts luminosity densities in 
to SFR densities. The resulting SFR densities are 
given in Table 3 and in Figure 3 are compared 
with others in the literature from both broad-band 
UV estimates and narrow line estimates (for the 
latter we go back to the original data and con- 
sistently use our conversions/extinctions for line 
luminosity densities to SFR densities). In gen- 
eral our results fall within the scatter of previous 
measurements given the large error bars from the 
small number statistics of our sample. Our Ha, 
H/3 points agree with the range of Balmer line es- 
timates at z < 0.5 from the literature, at z — 0.9 
our [Oil] estimate is a little low (even given the 
large one object error bars) especially compared to 
other [Oil] estimates. It is possible that our mea- 
sure is underestimated due to an incorrect choice 
of [Oil] /Ha ratio, for example if we adopt the 'un- 
biassed' [Oil] /SFR conversion of Rosa-Gonzalez, 
Terlevich, & Terlevich (2002) we can raise our 
[Oil] point by 0.35 dex. It is also possible that 
our assumed extinction value is incorrect. The ef- 
fect of increasing the extinction, for example to 
Ay = 2.0 mags would be greatest on our H/3 SFR 
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Table 3: Values calculated in determining the star formation rate density. Luminosities refer to equivalent 
dust-corrected Ha values derived from each actual line. 



2 


Line 


Log 
EiLi/V 
(erg sec 1 Mpc 3 ) 


Log 

<t>* 
(Mpc" 3 ) 


Log 
L* 

(erg sec -1 ) 


Log 
PL (TOTAL) 
(erg sec" 1 Mpc -3 ) 


Log 

PSFR 
(M© yr" 1 Mpc -3 ) 


Luminosity 
Function 
Limit 


0.38 


Ha 


39.36 


-2.55 


41.82 


39.41 


-1 70+ ' 14 

L - ' u -0.21 


.03 L* 


0.46 


H/3 


39.94 


-2.23 


42.15 


40.07 


-1 04+ ' 17 


.11 L* 


0.90 


[Oil] 


39.30 


-2.91 


42.30 


39.52 


i co+0.34 
1 -°°-0.30 


.23 L* 



density raising the estimate by 0.5 dex while the 
other estimates would be raised by 0.3 dex. This 
would make [Oil] agree better but then the H/3 
point would be highly discrepant. Similarly the ef- 
fect of assuming zero extinction would be to lower 
the points by the same amounts. Clearly then our 
SFR density measurements would still be incon- 
sistent: [Oil] would be highly discrepant. Inter- 
nal consistency gives us good reason to think our 
choices of [Oil] /Ha line ratio and extinction are 
reasonable, though the allowed range is large. 

The formal error bars also do not include the 
effect of large scale structure which are signifi- 
cant. Our small volume is equivalent to a sphere 
of radis 5-6 Mpc which is ~ the local blue galaxy 
correlation length. (Madgwick et al. (2003) finds 
r = 5.2-5.9 Mpc for star-forming galaxies.) For 
the volumes (~ 500-1000 Mpc 3 ) and number den- 
sities (0.002-0.02 Mpc~ 3 ) considered here we ex- 
pect fluctuations of 100-130% relative to the mean 
using the results of Somerville et al. (2004). The 
deviations of the points from the line in Figure 3 
are of this order and so large scale structure would 
represent an equally reasonable explanation of our 
scatter. 

5. Conclusions 

This work demonstrates we can select galaxies 
purely by line emission in a three-dimensional sur- 
vey, producing a pure narrow-band sample limited 
by line luminosity to constrain the cosmological 
SFR. One can see from Table 3 that our observ- 
able luminosity limit for Ha is almost two orders 
of magnitude lower than L* , indicating a small 
completeness correction. The luminosity densities 
calculated from this data are therefore very com- 
plete and provide a strong bound on the cosmo- 



logical total. 

Our study also finds no evidence for any new 
population of low-luminosity objects with strong 
emission lines which could have caused previous 
estimates of the SFR density to be seriously un- 
derestimated. Searching for such unexpected pop- 
ulations was a primary goal of our observations, 
we are reporting a null result for this search. 

Rather we find a SFR density history consis- 
tent (albeit with large errors from small num- 
ber/volume statistics) with earlier estimates from 
broad-band selected samples. More specifically 
our measured SFR density agrees with the pre- 
vious emission line determinations (both broad- 
band and narrow-band selected) but are higher 
than previous ultraviolet determinations which is 
in accord with earlier findings. Our low sub-L* 
line luminosity limits corresponds to seeing ~ 80- 
90% of the integrated SFR density at z - 0.4; 
even at z = 1 our flux limit corresponds to ~ 40- 
60% of the total star-forming light. This indicates 
that previous completeness estimates to broad- 
band surveys were correct and that we now have a 
reasonable census of the total star-formation rate 
density in the Universe out to z = 1. 

The success of narrow-band tunable filter sur- 
veys (both this work and those of Jones and Bland- 
Hawthorn (2001) and Hippelein et al. (2003)) au- 
gurs well for the next generation of tunable fil- 
ter 3D instrumentation. Galaxies can be located 
in data cubes purely from their line emission flux 
(i.e. not equivalent width) and recovered aperture- 
extracted spectra of general objects agree well 
with those taken using classical dispersive spec- 
trographs. For line emitting objects modest expo- 
sures on 4m telescopes allow one to survey sub- 
stantially below L* in volume-limited samples at 
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Fig. 2. — Equivalent Ha luminosity (dust- 
corrected) versus rcdshift of individual galaxies. 
Also shown are lower limits calculated from our 
flux limits. The solid line shows for reference the 
shape of the cosmological SFR density- z curve as 
parameterized by Cole et al. (2001) (extinction 
corrected) re-normalized to an Ha luminosity of 
10 4i.86 ergs s -i at z = o (this i s t he value of L* 

quoted by Gallego et al. (1995)). The dashed, hor- 
izontal line shows the equivalent luminosity cor- 
responding to SFR = const. = 1 A/ yr. _1 (We 
assume no redshift dependence of the Kennicutt 
L(Ha)/ SFR conversion.) 

high redshift. Our survey suffers from a limited 
volume surveyed using an instrument with a 5 
arcmin FOV; future wider field instrumentation 
should allow much bigger surveys. In the near- 
future we also expect such instrumentation tech- 
niques to be extended beyond the CCD cutoff to 
the near-IR and especially the J-band. Here the 
gain by observing in narrow wavelength slices is 
even greater than in the optical because the air- 
glow lines are even stronger. Here tunable filter 
techniques (Bland-Hawthorn & Nulsen 2004) on 
8m class telescopes would allow astronomers to 
find Ly-a emitters at z > 7 for even very mod- 
est Ly-a / SFR assumptions and thus allow us to 
experimentally probe the very first galaxies and 
the epoch of re-ionization. 
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z data of Cole et al. (2001) (extinction corrected 
version) . 
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